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Measurements of the production of forward jets from polarized proton collisions at ^/a — 500 GeV 
are reported. Our measured jet cross section is consistent with hard scattering expectations. Our 
measured analyzing power for forward jet production is small and positive, consistent with the 
existence of partonic orbital angular momentum. 
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The proton is a building bloek of matter, which is it- 
self built from elementary quarks (g) and gluons (g). Our 
understanding of the structure of the proton has become 
increasingly sophisticated since the advent of Quantum 
Chromodynamics (QCD), and a reason for this has been 
the quest to understand how the proton gets its intrin- 
sic spin from its constituents. The present view is that 
q or g orbital angular momentum (0AM) makes impor- 
tant contributions to the proton spin. Early indications 
of this came from large analyzing powers (Ajv) measured 
in the production of charged and neutral pions in colli- 
sions of transversely polarized protons at center-of-mass 
energy y/s = 20 GeV [l| . A^ is the amplitude of the spin- 
correlated azimuthal modulation of the produced parti- 
cles. A large A^ is not expected for pions produced with 
sufficient transverse momentum (pt) in coUinear pertur- 
bative QCD (pQCD) at leading twist, due to the chiral 
properties of the theory [2] . Measurements of large An 
for TT production at large Feynman-x {xp = '^Pz/Vs) 
prompted theorists to introduce spin-correlated trans- 
verse momentum (kx) in either the initial state (Sivers ef- 
fect jsl) or the final state (Collins effect [4|). For inclusive 
TT production, these effects cannot be disentangled. In 
contrast, in measurements of jets, defined as a collimated 
multiplicity of energetic baryons and mesons that are pro- 
duced in high-energy collisions, contributions from frag- 
mentation are absent and hence information about the 
scattered q or g can be inferred directly. In particular. 



An for jet production, direct photons or Drell-Yan pro- 
cesses is expected to arise only from the Sivers effect. 
The initial-state spin-correlated fcy is related by models 
to g and g 0AM. 

Large An for tt production at large xp in p^ + p col- 
lisions at -y/s < 20 GeV [l|, Iq have challenged under- 
standing, because cross sections are much larger than 
pQCD expectations, resulting in skepticism that tt pro- 
duction in these kinematics is from hard-scattering pro- 
cesses. Theoretical interest in understanding An for 
TT production has been revived by recent measurements 
[Tl |8| at ^/s > 62 GeV, where cross sections [9| are in 
agreement with pQCD. Furthermore, measurements of 
An for TT production in p^ +p collisions at -^/s > 62 GeV 
have been concurrent with measurements of transverse 
single-spin asymmetries (SSA) in semi-inclusive deep in- 
elastic scattering (SIDIS) [10| , where an electron or muon 
is inelastically scattered from a proton, whose spin is 
transverse to the lepton beam. Meson fragments of the 
struck q are found to have spin-correlated azimuthal 
modulations, whose amplitudes are understood by the 
Sivers and Collins effects, introduced to explain An for 
p^+p -^ TT + X. An alternative and complementary theo- 
retical approach based on collinear factorization ll| pre- 



dicts An involving twist-3 multi-parton correlations [1^ , 
and is expected to be related to the Sivers and Collins 
functions via transverse-momentum moments. However, 
an attempt at linking the different approaches using data 



from SIDIS and p^+p -^ n + X yielded a mismatch in the 
sign [l3|. Thus a consistent understanding of all trans- 
verse SSA in hard scattering processes is not yet within 
our grasp, but would greatly benefit from measurements 
of An for jet production in p^ + p collisions, since it re- 
ceives no contributions from the Collins effect. 

In this Letter, we report first measurements of cross 
sections and An for forward jet production inp^ +p col- 
lisions at -y/s=500 GeV. The measurement was conducted 
with the AivDY detector at the 2 o'clock interaction re- 
gion (IP2) of the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory. The primary de- 
tector components were two mirror-symmetric hadron 
calorimeter (HCal) modules that were mounted to face 
one of the beams at IP2 for the 2011 and 2012 RHIC runs, 
and spanned the pseudorapidity interval 2.4 < 77 < 4.0 
that is well shielded from single beam backgrounds by 
the cryostats of the ring magnets. 

Each HCal consisted of a 12-row x 9-column matrix 
of (10 cm)^ X 117 cm lead cells, each with an embed- 
ded 47 X 47 matrix of scintillating fibers [1^. For the 
2012 run, two 5-row x 2-column arrays were deployed 
above and below the beams to create an annular HCal 
with a central 20 x 20 cm^ hole for the beams. The HCal 
positions were measured by survey. The A^rDY appa- 
ratus also had a pair of 16-element scintillator annuli 
mounted symmetrically about IP2 to serve as a beam- 
beam counter (BBC) [15[, a pair of 7- row x 7- column lead 
glass detector arrays serving as small electromagnetic 
calorimeters (ECal) at a fixed (variable) position for the 

2011 (2012) run, a scintillator preshower array and a pair 
of zero-degree calorimeter (ZDC) modules [17| that faced 
each beam. A GEANT [l^ model of A^vDY was created, 
as described in Ref. [18[, and uses inputs from PYTHIA 
6.222 [19[, hereafter referred to as full simulation. 

p^ +p collisions were initiated at IP 2 at systematically 
different times in stores during the 2011 run to assess 
the impact of a third collision point on operations. The 
colliding beam luminosity at IP2 was measured by the 
vernier scan technique and resulted in ct=0.94±0.08 mb 
for the effective cross section of coincidences between the 
ZDC modules which were used to continuously monitor 
the luminosity. Results reported here were from 6.5 pb^^ 
of integrated luminosity during the 2011 run at ^/s—500 
GcV and 2.5 pb~^ of integrated luminosity during the 

2012 run at -^i = 510 GeV. The polarization of each 
beam was measured by a relative polarimcter at several 
times in each fill. The relative polarimeter was calibrated 
from measurements from an absolute polarimeter, result- 
ing in 0.526±0.027 (0.511±0.028) ^ for the average 
beam polarization for our jets at xp > {xp < 0) in 
the 2011 run. 

The data by AatDY is from 32-channel 70 MHz flash 
analog-to-digital (ADC) converters with 0.25 pC/count 
sensitivity and noise levels < 0.25 pC. Online pedestal 
corrections were made on board by field-programmable 
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FIG. 1: Tower multiplicity distributions for forward jets for 
data compared to full simulations for (left) jets from ^5=500 
GeV collisions, as used for the jet analyzing power; and for 
(middle) jets from v^=510 GeV collisions, as used for the jet 
cross section. (Right) Multiplicity of particles produced by 
PYTHIA 6.222 ^ that gives rise to the forward jet. 



gate arrays (FPGA). Pedestal-corrected ADC counts 
were then analyzed by the FPGA to derive an event trig- 
ger. The majority of the data were from our jet trigger 
that summed the ADC response from each HCal mod- 
ule, excluding the outer two perimeters of cells. This 
trigger is sensitive to electromagnetic (EM) and charged 
hadron fragments of jets. Events were also acquired from 
a minimum-bias trigger that required minimum charge 
from any clement of the annular BBC that faced each 
beam, sometimes with a collision vertex requirement, 
where the vertex is reconstructed by the FPGA from 
time-difference measurements. Some events were ac- 
quired when either ZDC crossed threshold as a way to 
tune bunch-crossing scalers used to monitor luminosity 
and the polarization of colliding beams. 

The HCal had individual cell gains adjusted prior 
to colliding beam operation based on their cosmic-ray 
muon responses. The first step in the offline analysis 
was to determine the absolute energy scale of the HCal 
modules by reconstruction of tt'' — >■ 77 from pairs of 
single-cell clusters, that had neighboring cells with en- 
ergy E' < 0.11 GeV, where E' is the incident photon 
equivalent energy. The tt" identification was confirmed 
by associating the singlc-ccU clusters reconstructed from 
our full simulations with primary photons generated by 
PYTHIA [l8|. Offline analysis also refined the relative 
calibration of all cells by a combination of 7r° reconstruc- 
tion and the matching of energy-deposition distributions 
from single cells between data and full simulation. For 
the jet analyses described below, the energy calibration 
was adjusted to account for the average difference be- 
tween hadronic and EM showers from full simulations. 
We used E = 1.12x£" — 0.1 GeV, where E is the equiva- 
lent incident energy measured by individual cells as used 
in the jet finders and E' is from tt*^ calibration. 
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FIG. 2: (left) Event averaged jet shape, corresponding to how 
the energy depends on R, the distance a tower is from the 
thrust axis in (ry, 0) space, (right) Correlation between jet xf 
and pt ■ 



Our final jet results use the anti-fcy jet algorithm (2l| 
with a cone radius of Rjet=0-7 radians in (77, 0) space, al- 
though we have also considered cone algorithms [22] ■ For 
each event, the iji, <J3i of each cell with E > Ethr [Ethr = 
0.25 GeV is used to make towers) is reconstructed from 
its surveyed position and the z position of the collision 
vertex for the event, calculated from the measured time 
difference between the two BBC annuli. The anti-fcy 
algorithm reconstructs the jet by a pair-wise merging 
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of tower clusters separated by di. 

(Rfj/R'j^f), when dy < 1/k'^-ior any i. Each tower has a 
transverse momentum fcr.i = Ei/cosh.{r]i), assuming zero 
mass for the incident particle. Pairs of tower clusters are 
separated by Rij = V('?* ~ ''h)^ + (<^« ~ '?^i)^- The merg- 
ing procedure is repeated until all towers are accounted 
for. A valid jet is considered to have \rjjet — ?7o| < drj 
and \4>jet — 0o| < d(j>, where 'rij(.t{4>jet) are computed from 
the energy-weighted averages of towers included in the 
jet within the acceptance centered at tjq, (J)q of half- width 
dr?, d<j). For the 2011 data, both ECal and HCal cells 
were considered. For the 2012 data with ECal positioned 
beyond the HCal acceptance, only HCal cells are con- 
sidered. Our cross section and analyzing power results 
are reported using 770=3. 25, d77=0.25 and d(j)=0.5. The 
HCal module to the left (right) of the oncoming beam 
has (/>o = (tt) . 

The tower multiplicity distributions for valid jets are 
shown in Fig. [1] This figure compares jets reconstructed 
from data to jets reconstructed from our full simulations 
oip+p collisions, where the GEANT response uses the in- 
dividual cell calibrations to produce simulated ADC val- 
ues, and the jet trigger is emulated by the same algorithm 
used by the FPGA for our measurements. In general, 
the simulation gives an excellent description of the data, 
consistent with minimal contributions from single-beam 
backgrounds, as determined from direct measurement, or 
from other unknown sources of energy deposition (under- 



lying event). Small increases in the HCal multiplicity for 

2011 are attributed to ancilliary material (e.g., cables, 
etc.) not included in GEANT, prompting us to use the 

2012 data for the jet cross section. The reconstructed jets 
have a broad tower multiplicity distribution whose mean 
value increases as Ej^t increases. Given that data and 
full simulation agree, we can then infer the distribution 
of particles in the jet by applying the anti-fc^ jet finder to 
detectable particles produced by PYTHIA. These parti- 
cle jets have similar multiplicity to those reconstructed in 
e"'"e~ collisions at ^/s « 10 GeV [23| and in fixed-target 
hadroproduction experiments [25| . Such low multiplicity 
jets are generally not accessible in hadron colliders be- 
cause their pT is too low. Forward detection results in 
large xp, so makes their detection possible. 

The towers included in the jets have their energy dis- 
tributed relative to the thrust axis in a manner that is 
typical of a jet (Fig. ^. Most of the energy is concen- 
trated near the thrust axis. As towers become increas- 
ingly distant from the thrust axis, on average they con- 
tribute little to the energy of the jet. The data is well 
described by our full simulation, although there are some 
indications that jets produced by PYTHIA 6.222 ^ 
have energy concentrated closer to the thrust axis rel- 
ative to our measurements. Also shown in Fig. [2] is the 
correlation between xp and pt for our jet events. The 
di] requirement strongly correlates these two kinematic 
variables. 

The jet energy scale was established by comparing 
tower jets reconstructed from the full simulation to par- 
ticle jets reconstructed from PYTHIA, and resulted in 
the hadronic compensation described earlier. A check of 
the energy scale was made for 3-jet events in the data 
by the centroid of a peak attributed to T(1S)— >■ 3g [23J . 
This peak has statistical significance of 3.5a. The cen- 
troid of the peak depends smoothly on Rjet used in the 
anti-fcy algorithm, as do our measures of the jet-energy 
scale from simulation. 

The forward jet production cross section was measured 
by scaling the number of reconstructed jets by the mea- 
sured integrated luminosity and correction factors de- 
scribed here. For jet triggers, there is a trigger efficiency 
i^trig) dominated by the variation of ry along the colli- 
sion vertex distribution, etrig is evaluated as a function 
of jet energy from the full simulation, and is checked by 
comparing invariant jet cross sections from jet-triggered 
events to cross sections determined from the minimum- 
bias trigger. The jet detection efficiency (ejet) is deter- 
mined from the ratio of number of tower jets within the 
acceptance to the number of particle jets within the ac- 
ceptance, ejet is checked by systematically varying the 
acceptance and assessing the stability of the resulting in- 
variant cross section. Sources of systematic uncertainty 
are (a) values oi etrig and edet', (b) time-dependent effects 
from cither HCal gain stability or from beam conditions; 
(c) jct-encrgy scale uncertainties; (d) luminosity normal- 
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FIG. 3: Invariant forward jet cross section measured at 
rj — 3.25. Our measurements are compared to predictions 
by PYTHIA. 



ization; and (c) jct-findcr parameters (Rjet, Ethr) that 
also probe underlying event contributions. The result- 
ing distribution of forward jet cross section as a func- 
tion of energy is shown in Fig. |3l Our results are com- 
pared to PYTHIA 6.222 ^ and 6.425 ^ predictions 
for anti-ZcT jets reconstructed from stable particles that 
are within the detector acceptance. In Ref. j22l | it was 
shown that PYTHIA predicts forward jets are from par- 
tonic hard scattering. PYTHIA 6.222 is prior to tunings 
based on Tevatron data which resulted in later versions 
(e.g., 6.425) used by the LHC. Versions of PYTHIA that 
predate tunings for the LHC are known to accurately de- 
scribe large xp tt^ production [26| , and are known to lose 
accuracy for more complicated multi-particle correlations 

The forward jet An is measured by the cross-ratio 
method, that cancels through second order luminosity 
and detector asymmetries by combining spin up (f) 
and spin down (J,) measurements from nominally mir- 
ror symmetric beam-left (L) and beam-right (i?) hadron 
calorimeter modules. The jet analyzing power is 
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Each fill has a pattern of spin directions for bunches 
of beam injected into RHIC. A specific crossing of 
bunches from the two rings is the remainder after di- 
viding the RHIC clock count for an event by 120. The 
bunch-crossing distribution has characteristic holes that 
correspond to missing bunches from one or the other 
beam. The pattern of polarization directions for that fill 
recorded at AjvDY originating from information broad- 
cast by RHIC, is then used to accumulate NJ^Z^. in the 
analysis. Since the RHIC broadcast information speci- 
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FIG. 4: Analyzing power for forward jet production. Jets 
are reconstructed with the anti-fcr algorithm using Rjet = 
0.7. Preliminary results Q^ reported comparable An with 
the mid-point cone algorithm. Systematic error estimates are 
described in the text, and do not include scale uncertainty 
from the beam polarization measurements. 



fies polarization directions at the polarized ion source, 
we rely on the measurement of spin asymmetries for far- 
forward neutron production measured by the ZDC, where 
the Ajv was previously measured 
An is measured with the proper sign. 



28[, to ensure the jet 



Our measured forward jet An is shown in Fig. 3] One 
check for systematic effects was to fit the spin asymmetry 
(e = PbeamAN) mcasurcd in each jet (xp) bin for each 
RHIC fill by a constant. The resulting x^ per degree of 
freedom from these fits is close to unity, and is consistent 
with the statistical errors, meaning the systematic errors 
are small. A more quantitative check for systematic ef- 
fects was to establish if an effectively unpolarized sample 
oi p + p collisions had An consistent with zero. This was 
accomplished by a random reversal of the spin direction 
for half of the filled bunch crossings. The mean value of 
e for ^100 random spin direction patterns had values 
10^^ < e < 10~* resulting in the systematic uncertainty 
estimate of 2 x 10~* for the jet An- The systematic er- 
ror estimates in Fig. |4] are estimated by varying the jet 
finder and valid jet parameters. Our jet An measurement 
is limited by statistics. Our measured small and positive 
jet An is naively expected because ^Ar(7r"'")sa — AAr(7r~), 
thus giving cancelling contributions from tt* in a jet. 

In conclusion, we have made first measurements of for- 
ward jet production in p^ +p collisions at -^5=500 GeV. 
Our measured cross section is consistent with dominant 
contributions from partonic hard scattering, even though 
the transverse momentum for the produced jets is small 
(2 < pt < 10 GeV/c). We have measured the analyz- 
ing power for forward jet production, and find it to be 
small and positive. Our measurements constrain knowl- 
edge [29| of Sivers functions, that are related to parton 



0AM through models. It remains the case that the most 
definitive experiment to test present understanding is a 
measurement of the analyzing power for DY production. 
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